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Abstract 



We present theoretical analyses and numerical simulations for the adhesion-induced 
phase separation of mult i- component membranes with two types of ligand-receptor 
complexes (junctions). We show that after integrating all possible distributions of 
the junctions, the system can be regarded as a membrane under an effective external 
potential. Mean field theory and Gaussian approximation are used to analyze the 
effective membrane potential and we find (i) The height difference of the junctions 
is the main factor that drives phase separation at sufficiently large junction height 
difference, (ii) In the two phase region far from the mean- field critical point, because 
of the higher entropy associated with the softer junctions, phase coexistence occurs 
when the effective binding energy of the more rigid junctions is higher, (iii) In the 
two phase region near the mean- field critical point, the shape of the effective potential 
shows that the phase coexistence occurs when the effective binding energy of softer 
junctions is higher. The effect of junction density on the critical point is studied by 
Monte Carlo simulations, and the result shows that phase separation occurs at larger 
junction height difference as junction density of the system decreases. 
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Notation 

Material parameters 



symbol 


dimensionless form 


physical meaning 


K 




bending rigidity 


7 


la a? 
K 


surface tension 


X a 


X a a 2 

K 


elastic constant of type-a junction 


K 


P 

v a 


natural length of type-a junction 


Ebcx 


Eboc 


binding energy of type-a junction 




4>a 


density of type-a junction 


/J, a 


P-a 


chemical potential of type-a junction 


a 




lattice constant in the x — y plane 


h 




unit length in the ^-direction 


k B T 




energy unit 



V 



Defined parameters 

£„ = ^ 
A h = ^ 
A W 

a K 

A+ = Ai + A 2 
A_ = Ai - A 2 

A - A - 

E eff+ = E efh + E efh 
E eff- = E efh - E eff2 

9 = A h 2 A + 



the average height of the two junctions 
length difference of the two junctions 
dimensionless junction elastic constant 
sum of the junction elastic constant 
difference of the junction elastic constant 
junction flexibility difference 
effective hinging energy 
sum of the effective binding energy 
difference of the effective binding energy 
junction height difference 
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Chapter 1 
Introduction 



Membrane adhesion has been studied theoretically and experimentally [1]- [1 1] . It 
plays an important role in many biological processes ranging from embryological de- 
velopment, repair of tissue, immune response [12], to pathology of tumors [IS]- In 
general, the adhesion is mediated by some kinds of lock-and-key molecular complexes 
(for simplicity, they are called junctions). An example of membrane adhesion is the 
formation of immunological synapse between a T lymphocyte (T cell) and an antigen- 
present cell (APC), a key event of the immune response ^2]- in this example a highly 
organized pattern in membrane adhesion region where the TCR/MHC-peptide com- 
plexes aggregate in the center with a peripheral ring composed of the LFA-l/ICAM-1 
complexes is observed. 

Qi et. al. [14] studied the immunological synapse pattern formation with a 
theoretical model based on a Landau-Ginzburg free energy and a set of coarse- 
grained react ion- diffusion equations and suggested that the dynamics of immuno- 
logical synapse pattern formation is a spontaneous self-assembly process. Later, Ray- 



1 



chaudhuri et. al. PH] developed an effective membrane model which is derived from 
the reaction-diffusion equations in ^3] and the condition of forming immunological 
synapse is studied by the mean field, Gaussian, and renormalization group theories. 

In [Hi] Lipowsky et. al. studied membrane adhesion in the presence of one type 
of junctions and one type of repellers (repulsive molecules) in their Monte Carlo Sim- 
ulation. The results show that the phase separation depends on the height difference 
and the concentration of two types of molecules. Later in jTZJ , they considered mem- 
brane adhesion in the presence of two different types of junctions which resembles 
T cell-APC adhesion. Their study indicated that the phase separation is driven by 
the height difference between two types of junctions, but the formation of target- 
patterned immunological synapse has to be assisted by the motion of cytoskeleton. 

We study the general case of two membranes binding to each other due to the pres- 
ence of two types of junctions. Similar system has been studied in [TH], where Chen 
considered the effect of junction flexibility difference on phase separation and devel- 
oped an equilibrium statistical mechanical analysis which provided a phase diagram 
in mean field level. In this thesis, we re-examine the system studied in jTHJ with an 
effective membrane model that is closely related to ^5] • Our goal is to provide a com- 
plete picture for the phase separation induced by membrane adhesion which includes 
the effect of junction height difference, junction flexibility difference, and thermally 
activated membrane height fluctuations. This thesis is organized as follows. Chapter 
2 introduces the model and the effective potential approach, meanwhile mean field 
theory and Gaussian approximation are introduced to study the phase diagrams of 
the system. Chapter 3 discusses the simulation. Chapter 4 summarizes this thesis. 
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Chapter 2 
The Model 



In this chapter, we introduce the Hamiltonian of the system under consideration 
and derive the effective potential acting on the membrane. Through the analysis of 
effective potential in the mean-field approach and Gaussian approximation, a lot of 
the physics of this system are revealed. 

The main compositions of biomembranes are amphiphilic lipids with hydrophilic 
heads and hydrophobic tails. Besides the lipids, there are many kinds of proteins 
and carbohydrates. For simplicity, in our model we consider membranes containing 
only lipids, ligands, and receptors. The adhesion between a substrate-supported 
membrane and a membrane floating in the solvent is induced by the formation of 
two types of specific ligand-receptor complexes (in the rest of this thesis, we shall 
refer them as "junctions"). The geometry of the system under consideration is shown 
schematically in Fig. 12.11 The height of the upper membrane at r is denoted as h(r), 
where r = (x, y) is a two-dimensional planar vector. The junctions which mediate the 
membrane adhesion are formed by non-covalent adhesion between type-a (a is 1 or 2) 
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ligands and type-a receptors. The density of type-a junctions at r is denoted as $ Q (r), 
and the binding energy of a type-a junction is Es a - These ligand-receptor bonds are 
non-covalent bonds composed of several hydrogen or van der Waals- like interactions. 
Typical binding energy for a ligand-receptor pair is on the order of 10 ksT. Our 
model geometry directly applies to the experimental set up in jT^j. In Appendix 
A, we show that by suitable transformation, the system with two membranes in the 
solvent is equivalent to to the system described above. 



The energy of this system at any instance has several contributions: (i) the curva- 
ture elasticity and surface tension of the membrane, (ii) the binding energy of the 
ligand-receptor pairs, (iii) the deviation of the ligand-receptor pairs from their nat- 
ural conformation. Therefore the simplest coarse-grained Hamiltonian that includes 
all the above effects has the following form 



Here k is the bending rigidity of the membrane (typical value of k is 2 — 70 x 10 -20 N-m 
[19]). 7 is the surface tension of the membrane (typical value of 7 is 24 x 10~ 6 N/m [20j. 
In this simple model it is assumed that k and 7 are independent of the composition of 
the membranes, and the direct interactions between the junctions, free ligands, and 
free receptors are neglected. The nonspecific interactions between the membranes are 



2.1 The Hamiltonian 



H 
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Figure 2.1: The geometry of the system. The system is composed of a substrate- 
supported membrane in the bottom and an upper membrane floating in the solvent. 
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Table 2.1: Typical values of material parameters measured for resting cells 



parameter symbol typical value 



bending rigidity k 2 - 70 x 1(T 20 iV ■ m [TP] 

surface tension 7 24 x 1CT 6 N/m [2D] 

junction elastic constant A a 10~ 5 — 1 iV/m [21] 

junction length h a 10 — 30 nm 

junction binding energy E Ba ~ 10 /c#T 

number of or ligands (receptors) 10 5 /cell [2T| 



also neglected. When a type-a junction is at r, the interaction energy between the 
membranes acquires a minimum at h = h a (the natural height of a type-a junction, 
typical value 10 - 30 nm [2~Tj ). and the coupling term Y^ a =i ^a( r )lf (^( r ) ~ ha) 2 
comes from the Taylor expansion of h around h a , where X a is the flexibility of a 
type-a junction (typical value of A a is 10~ 5 — 1 N/m [21] )• The typical values of 
the material parameters mentioned above are listed in Table. 12.11 To proceed our 
discussion, we put the variables h(r) and $ a (r) on a two-dimensional lattice with 
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lattice constant a. The size of a is chosen to be the smallest length scale beyond 
which the continuum elastic description of membranes breaks down (a = 6 nm |22j). 
The energy unit is chosen to be fc^T, where the unit length in the x — y plane is a 



and ho = a/ y f^f m ^ ne z direction. That is, Ebu = Esa/ksT, l{ = h(r)/h and 
4> a (i) = a 2 $ a (r). So the non-dimensional form of the Hamiltonian is 



N 

2 re 

0a W 



i=l ^ 
2 

+E- 



(2.2) 



2 re 

Where the discretized Laplacian A^ in two dimension is given by A<^ = Ad£ x , y = 
£x+l,y £x—l,y £x,y+l -44,j/, and the discretized gradient V is = V£ x>y = 

^(i x +i,y—£x-i,y)x + ^(£ x ,y+i—£x,y-i)y- Since the calculation in the canonical ensemble 
(system with fixed total number of ligands and receptors) is cumbersome. Thus, for 
simplicity the analysis is carried out in the grand canonical ensemble and we introduce 
the chemical potential of type-a junction, \i a (/i a = ksTfla). Here we apply the hard- 
core repulsion between the junctions, consequently each site can be occupied by lipids 
or one type-1 junction or one type-2 junction. Therefore, the partition function of 
the system can be written as 

r N r 2 2 

Z = V[£] e~ H ^ ft 1 + e^'H^fc-*) 3 + e ^// a -^(A-fe) a ( 2 . 3 ) 
V[£] e-G^+^/z), (2.4) 



where -E e // a = E Ba + fx a is the effective binging energy which represents the combi- 
nations of binding energy of ligand-receptor pairs and entropy lost of free ligands and 
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free receptors and 

TV 



i=l 



l 7 a 2 



2(W+2^-(VA) 



(2.5) 



is the elastic energy of the membrane and 



TV 



^e// = -E ln 



i=l 



a=l 



(2.6) 



is the effective potential acting on the membrane by summing over all possible dis- 
tributions of junctions. V e ff(i) is the effective potential at site i. The total number 
of type-a junctions N a in the system can be expressed as the derivative of the Grand 
potential with respect to the chemical potential. 



G = -k B T\nZ, N a = 



dG 

dn a 



E<- 



iV 2 



(2.7) 



Since typical experiments (either living cells or artificial membranes) are carried out 
for systems with fixed total number of ligands and receptors. At the end of our anal- 
ysis we connect E e ff a to the corresponding junction binding energies and densities of 
ligands and receptors in this system, and compare these values with typical biological 



systems . 



2.2 Phase Diagram: Zero Fluctuation 

As has been shown in the previous section, the equilibrium membrane height confor- 
mation of the system is the same as that of a membrane with effective Hamiltonian 



H eff = H e£ + U e ff 



(2.8) 



N 

E 



Table 2.2: Dimensionless parameters introduced in Sec. 2.2 



— tl+t-2 



A h = ^ 



A + = Ai + A 2 , A_ = Ai - A 2 

E eff+ = Eeffi + E e ff 2 , E e ff- = E effl - E e ff 2 

9 = A ft 2 A+ 



For convenience, several dimensionless parameters are defined in Table. 12.21 Sub- 
stitute these parameters into Eq. (|2.6|) . V e ff(i) can be expressed as a function of 



= z, 



E e.f.f+-9/2 

1 + e a A(z) 



(2.9) 



where 



A(z) 



, Ee// V" V2 -<^[('+i) 2 -i] + e - fie// - 2 ' gV2 -^[(^-D 2 -i] 



(2.10) 
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The equilibrium membrane height is determined by minimizing V e ff(i) with respect 
to z, i.e., 

E e f f+ -g/2 

dV<ff(f) _ e » dA(z) _ 

From Eq. (|2.11j) . phase separation occurs when dVel d l^ = has three roots (one 
maximum and two degenerate minimums). Evidently, the value of z that minimizes 
V e ff(i) is independent of E e ff + . Therefore, the phase boundary is also independent 
of E e ff + , and only depends on A, g, and -E e //-- However, it will become clear in the 
next section that this only holds in the mean field theory. 

2.2.1 Symmetric Case 

In the symmetric case (A = 0, i.e., both types of junctions have the same flexibility), 
the effective potential becomes 

= - In L^(e^iK^+e-^iN 2 -'])| . (2.12) 

In this case, the minimum of V e ff(i) depends on g and E e ff_. Since when £^ e //- — 0) 

A(z) = e-^ +1 M + e"^- 1 ) 2 - 1 ] (2.13) 

is an even function of z. Phase separation occurs at sufficiently large g, and there 
are three solutions for dA ^ = 0. For small g, there is only a triple root at z = 
which corresponds to £ = +£2)- And we find that the critical point is located at 
(E e ff- = 0, g = 2). The shape of V e ff(i) for two typical cases when phase separation 
occurs are shown in Fig. 12.21 The parameters for both lines are Ai = A2 = 0.2, 
E effl = E eff2 = Mi = 15, for solid line £ 2 = 30 (g = 22.5), for dashed line £ 2 = 22 
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Figure 2.2: The shape of V e ff(i) (effective potential) as a function of £ for two typical 
cases when phase separation occurs. The parameters for both lines are Ai = A 2 = 0.2, 
E effl = E eff2 = 1, £i = 15, for solid line £ 2 = 30 (g = 22.5), for dashed line £ 2 = 22 
(g = 4.9). The barrier height between the two minimums of V e ff(i) is denoted as A 
which is non-zero when phase separation occurs. The minimums are near £\ and £ 2 , 
the maximum is near £ = (£ x + ^)/2. It also shows that A becomes smaller when 
£■2 — £\ decreases. 
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Figure 2.3: The A — g diagram of symmetric case. For different -E , e //+ (by setting 
E e ffi = E e ff 2 ), the location of the critical point is g = 2 which is independent of 
E e ff + . And A decreases as -E , e //+ decreases. 
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(g = 4.9). The barrier height between the two minimums of V e ff(i) is denoted as A 
which is non-zero when phase separation occurs and becomes smaller when £2 — i\ 
decreases. Numerical calculation leads to the critical point locates at (E e ff_ = 0, 
g = 2). This result agrees what we obtain previously. Furthermore, Fig. 12.31 shows 
that for different E e ff + (by setting E e ffi = E e ff 2 ), the critical value of g is always 2. 
It confirms that the critical point in the mean field theory is independent of E e ff + . 
And we find that the values of A decrease as E e ff + decreases. That is, the depth of 
the potential minimums becomes shallower for smaller effective binding energy. From 
the above analysis, we know that the phase separation is driven by junction height 
difference, and this separation only occurs when g is sufficiently large. 

2.2.2 Asymmetric Case 

Let us consider the more general case A ^ 0, i.e., different types of junctions have 
different flexibilities. In this case, the symmetry in A = case no longer exists. 
Therefore we expect that the phase coexistence curve is no longer located at E e ff_ = 
line. An analytical expression of the phase coexistence curve can be found for 
I A I <gC 1 by expanding dA ^ around A = 0, where A(z) is taken from Eq. ()2.10|) . The 
critical point for A 7^ locates at the triple root of dA(z)/dz = 0. A straightforward 
calculation similar to ^H] leads to the position of the critical points for small A, 



9 = 2(1 



4 



(2.14) 



E, 



eff- - ~ 



A + C(A 3 ). 



(2.15) 
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The critical value of g decreases as the junction flexibility difference increases, and 
E e ff_, A have opposite sign. We present the phase coexistence curve in two differ- 
ent ways. Fig. 12.41 shows phase boundaries for E e ff + = 2, A + = 0.3, and A = 1/3 
(square), 4/15 (circle), and 1/6 (triangle). Fig. 12.51 shows phase boundaries for 
E e ff+ = 2, Ai = 0.2, A+ = 0.3, and A = 1/3 (square), 1/7 (triangleleft), and 1/15 
(triangleright). Both figures indicate that near the critical point, phase coexistence 
curves shift away from E e ff_ = line due to the junction flexibility difference. Fur- 
thermore, at large g, the phase boundaries are close to E e ff_ = line. To understand 
the physics of the phase boundary in small g region, we plot V e ff(i) as a function of £ 
in Fig. EIHlwith E effl = E eff2 = 1, Aj = 0.2, A 2 = 0.1 (A = 1/3), t x = 15, and £ 2 = 22 



(g = 3.657). The solid curve is K//, dotted curve for — In 



1 + e E effl -^-(e-h) 2 



which 



is the effective membrane potential without the contribution of type-2 junction, and 



dashed curve for — In 



1 + e E ef! 2-±t(t-t2? 



which is the effective membrane potential 



without the contribution of type-1 junction. At small g, the phase boundary shifts 
away from -E e //- = line because the minimum energy of the more rigid junctions 
is "lowered" by the softer junctions. Thus, in the region where the height difference 
is not very large, phase coexistence occurs when the effective binding energy of the 
more rigid junctions is smaller (the effective binding energy of the softer junctions is 
larger) . 
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Figure 2.4: Mean field phase boundaries for E e ff + = 2, A + = 0.3, and A = 1/3 
(square), 4/15 (circle), and 1/6 (triangle). We find that near the critical point, phase 
coexistence curves shift away from E e ff- = line due to the junction flexibility 
difference. Furthermore, at large g, the phase boundaries are close to E e ff_ = line. 
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Figure 2.5: Mean field phase boundaries for E e ff + — 2, A x = 0.2, A + = 0.3, and 
A = 1/3 (square), 1/7 (triangleleft), and 1/15 (triangleright). We find that near the 
critical point, phase coexistence curves shift away from = line due to the 

junction flexibility difference. Furthermore, at large g, the phase boundaries are close 
to E e ff- = line. 
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Figure 2.6: Effective potential as a function of £ in the small g region with -E e //i = 
E eff2 = 1, Ax = 0.2, A 2 = 0.1 (A = 1/3), l x = 15, and £ 2 = 22 (g = 3.657). The 



solid curve is V e ff, dotted curve for — In 



1 + eEeffl-^V-ll? 



which is the effective 



membrane potential without the contribution of type-2 junction, and dashed curve 



for — In 



1 + e^" 2 



which is the effective membrane potential without the 



contribution of type-1 junction. The softer junctions are longer than the more rigid 
junctions. The minimum energy of the more rigid junctions is "lowered" by the the 
softer junctions. 
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2.3 Gaussian Approximation 

In the mean field analysis, we assume that the membrane height is chosen to mini- 
mize V e ff, phase coexistence occurs when there are two degenerate minimums, and 
the critical point is located at where barrier height between the degenerate mini- 
mums vanishes. In this section, we introduce Gaussian approximation to study the 
fluctuation effects on phase boundaries. Since in equilibrium, the membrane height 
fluctuates around the potential minimums, the free energy for a membrane with aver- 
age height at £ m i na (a minimum of V e ff) in the Gaussian approximation is calculated 
by expanding V e // around £ m ma, i-e., 



- In J p[£] e ~[^+/ d2r ^//(^^ (2.16) 



k B T 

Integrating out these Gaussian fluctuations, the free energy of a'th minimum per unit 



area is 



F* 1 [ (27T) / 27T 

V e ff(£ mina ) - — / dq In, — + const. (2.17) 



k B T CJJ ^ nu »' 4tt J ^ y g 2 + Ag + 
After comparing the free energy per unit area of the two minimums of V e ff, we plot 
the phase coexistence curves for several choices of parameters in Fig. 12 .71 and Fig. 12.81 
In large g region, because the softer junctions can be stretched or compressed easier 
than the more rigid junctions and have higher entropy, comparing to mean field theory 
predictions, phase coexistence curve moves toward to higher E e ff\ (E e ff 2 ) when A > 
(A < 0). In Fig. 12.71 phase boundaries for several systems with E e ff+ = 2, A + = 0.3, 
and A = 2/3 (circle), 1/3 (square), and 1/6 (triangle) are shown. In large g region, 
because the softer the junctions are, the higher the entropy they have. Thus, the phase 
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boundary shift is more significant for systems with lager junction flexibility difference. 
Fig. 12.81 shows the phase boundaries for systems with A = 1/3, -E e //+ = 2, A + = 0.3, 
and E e ff + = —4 (lozenge), -2 (circle), 2 (triangledown) , and 10 (triangleleft). In the 
large g region, the phase coexistence curves are located farther from the E e ff_ = 
line (the phase boundary predicted by the mean field theory) for systems with smaller 
E e ff + . It is because the entropic effect of the softer junctions is smaller when the 
total junction density is higher (i.e., -E e //+ is larger). Fig. 12.71 and 12.81 also indicate 
that in small g region, because the location of critical points are not altered in the 
Gaussian theory, therefore the phase boundary goes to the mean field critical point 
as g decreases. 
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Figure 2.7: Phase boundaries in the Gaussian theory for E e ff + = 2, A + = 0.3, and 
A = 2/3 (circle), 1/3 (square), and 1/6 (triangle). Because the softer the junctions 
are, the higher the entropy they have, comparing to the more rigid junctions. The 
phase boundaries in the large g region shift away from -E e //- = line (mean field 
phase boundary), and the shift is more significant as |A| increases. The critical points 
in the figure are taken from the result of the mean field theory, because the location 
of critical points are not altered in Gaussian theory. 
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Figure 2.8: Phase boundaries in the Gaussian theory for A = 1/3, E e ff+ = 2, A + = 
0.3, and -E'e//+ = _ 4 (lozenge), -2 (circle), 2 (triangledown) , and 10 (triangleleft). 
The phase boundaries shift away from E e ff_ = line (mean field phase boundary), 
and the shift is more significant as -E , e //+ decreases. It is because the entropic effect 
of the softer junctions is smaller when the total junction density is higher (i.e., -E e //+ 
is larger). The critical points in the figure are taken from the result of the mean field 
theory, because the location of critical points are not altered in Gaussian theory. 
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2.4 Summary 



In the mean field analysis, the critical point and phase coexistence curve is indepen- 
dent of the value of E e ff + . For symmetric case (A = 0), the critical point is located at 
(E e fj- = 0, g = 2). The phase coexistence curve lies on the E e ff- = line starting 
from g = 2. In this case, the phase separation is driven by the height difference of 
the junctions. For asymmetric case(A ^ 0), the critical value of g decreases as the 
junction flexibility difference increases, and the phase boundaries shift away from the 
E e ff- = line in smaller junction height difference region due to the minimum of 
the more rigid junction is "lowered" by the softer junction. In large junction height 
difference region, the phase boundaries are close to the E e ff_ = line. 

In the Gaussian theory, we find that because the state associated with membrane 
height closer to the natural length of the softer junctions has higher entropy, as a 
result, comparing to the mean field predictions, in large g region the phase coexistence 
curves in the Gaussian theory move toward higher -E e //i (-E e //2)when A > (A < 0). 
The phase boundaries shift is more significant when the junction flexibility difference 
(|A|) increases or the total junction density (E e ff + ) decreases. In small g region, 
because the location of critical points are not altered in Gaussian theory, thus the 
phase boundary goes to the mean-field critical point as g decreases. 
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Chapter 3 



Numerical Simulation 



3.1 Monte Carlo Simulation 

Monte Carlo simulation is a general name for simulations which use random sequences. 
It is a way to perform statistical sampling experiments on a computer. It is named 
after the casino city "Monte Carlo" in the Monaco principality. 

3.1.1 Metropolis Algorithm 

Metropolis algorithm is applied to our simulations. In the algorithm, if the system is 
at state i at some instance, first one randomly chooses another state j as the possible 
new state of the system. Let Pi and Pj be the equilibrium probability for state i and 
j to occur at temperature T. The energy difference between state % and state j is 
AH = Ej — Ei. If AH < 0, then the system changes its state to state j. If AH > 0, 
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then the computer generates a random number 1Z G [0, 1]. If 

n<S = ^- = e- AH , (3.1) 

then the system changes its state to state j, otherwise the system keeps at state i. It 
can be proved that in the long time limit the simulation generates a series of states 
that obey Boltzmann distribution [23J. 

3.1.2 Monte Carlo Steps 

In our simulations, the energy of the system depends on the conformation of the 
membrane. In each MC step, we repeat the following procedure N = L 2 times. First, 
randomly choose a site {i={x, y)). The trial conformation of the membrane is the 
same as the current state except £ Xjy — > £' = £ XtV + A£, where — 5 < A£ < 5. This 
choice corresponds to a hight change of maximum size ~ 5 nm. If AH < 0, £ x , y is 
replaced by £' x . If AH > , then generate a random number 1Z and compare to 
e -AH £ Q decide whether the height of membrane should be changed or not. We also 
choose periodic boundary condition in all our simulations. 

3.1.3 Snapshots in MC Simulations 

The parameters in the simulations are listed in Table. 13.11 Two snapshots of the 
contours of the membrane height with -E e //i = -E e //2 = 1, A-i = A 2 = 0.2, £\ = 15, 
and £2 = 21 (i.e., g = 3.6) are shown in Fig. 13.11 and Fig. 13.21 The parameters 
are chosen such that the system is in two phase region, far from critical region. In 
Fig. 13.11 the initial membrane height is set to be £1, the natural length of type-1 
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Table 3.1: Parameters in simulations 



symbol values 



T 0.0025 
Ai, A 2 0.2 

e 1 is 



junction. The height of the black patches is smaller than £±, the height of the dark 
gray patches is in the interval between £\ and \{£\ + £•})■ In Fig. 13 .2[ the initial 
membrane height is set to be £2, the natural length of type-2 junction. The height 
of the light gray patches is in the interval between \[£\ + £2) and £2, the height of 
the white patches is larger than £2. Both figures show that in two phase region, 
the equilibrium membrane height depends on the initial condition (i.e., the natural 
length of type-a junction). It is because the barrier height between the potential 
minimums is so large comparing to the thermal fluctuation. Thus, the distribution of 
the equilibrium membrane height is near the initial membrane height. Fig. 13.31 and 
Fig. 13.41 are snapshots for -£7 e //i = E e jf 2 = 1, A x = A 2 = 0.2, £\ = 15, and £ 2 = 20 
(g = 2.5) on a 64 x 64 lattice. The set of parameters corresponds to a system near 
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Figure 3.1: Snapshot of the contours of the membrane height with -E e //i = E e ff 2 = 1, 
A 1 = A 2 = 0.2, l x = 15, and £ 2 = 21 (g = 3.6) on a 64 x 64 lattice. The parameters 
are chosen such that the system is in two phase region, far from critical region. The 
initial membrane height is set to be £ ± , the natural length of type-1 junction. The 
height of the black patches is smaller than £ 1; the height of the dark gray patches 
is in the interval between t\ and \{t\ + h)- That is, in this region the equilibrium 
membrane height is near the natural length of type-1 junction, the initial membrane 
height. 
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10 £0 SO 40 SO £0 



Figure 3.2: Snapshot of the contours of the membrane height with -E e //i = E e ff 2 = 1, 
Ai = A 2 = 0.2, i\ = 15, and i 2 — 21 (g — 3.6) on a 64 x 64 lattice. The parameters are 
chosen such that the system is in two phase region, far from critical region. The initial 
membrane height is set to be £ 2 , the natural length of type-2 junction. The height 
of the light gray patches is in the interval between + £ 2 ) and £ 2 , the height of 
the white patches is larger than £ 2 . That is, in this region the equilibrium membrane 
height is near the natural length of type-2 junction, the initial membrane height. 
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critical point. In Fig. the initial membrane height is set to be the natural 
length of type-1 junction. In Fig. 13.41 the initial membrane height is set to be £2, the 
natural length of type-2 junction. In both figures, the black, dark gray, light gray, and 
white color represent the membrane height £ correspond to £ < £1, £\ < £ < 7^1 + ^2), 
\[£\ + £2) < £ < £2, and £ > £2 respectively. In this region, the equilibrium membrane 
height is independent of the initial condition because the thermal fluctuations make 
the membrane height change from shorter (higher) junction phase to higher (shorter) 
junction phase. The above figures show the snapshots of the equilibrium membrane 
height for large g (far from the critical region) and small g (near the critical region). 
However, it is very difficult to locate the critical point by looking at the snapshots. 
Thus, we introduce a systematic method "Binder cumulant" to determine the critical 
points. 
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Figure 3.3: Snapshot of the contours of the membrane height with -E e //i = E e ff 2 = 1, 
Ai = A 2 = 0.2, i x = 15, and £ 2 = 20 (g = 2.5) on a 64 x 64 lattice. The parameters 
are chosen such that the system is in one phase region, near critical region. The initial 
membrane height is set to be £±, the natural length of type-1 junction. The black, 
dark gray, light gray, and white color represent the membrane height £ correspond to 
i < £ u i x < t < \{i x + £ 2 ), + £ 2 ) < t < e 2 , and £ > e 2 respectively. 
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Figure 3.4: Snapshot of the contours of the membrane height with -E e //i = ^e//2 = 1, 
Ai = A 2 = 0.2, i-y = 15, t 2 = 20 (g = 2.5) on a 64 x 64 lattice. The parameters are 
chosen such that the system is in one phase region, near critical region. The initial 
membrane height is set to be £2, the natural length of type-2 junction. The black, 
dark gray, light gray, and white color represent the membrane height £ corresponds 
to £ < £ u £ x < £ < \{t x + £ 2 ), \{t x + £ 2 ) < £ < £ 2 , and £ > £ 2 respectively. 
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3.2 The Binder Cumulant and Critical Point 

In last section, we have shown that it is difficult to determine the critical point 
quantitatively by looking at the snapshots of the simulations. Thus, we introduce the 
Binder cumulants C 2 and C4 which are defined as 

C ~ ( " 2) C - ( " 4) (3 2) 

where z — -h YuiLi z i ls ^ ne s P a tial average of the order parameter, Z{ = t- L — £0, and 
(O) is the thermal averages of O. For g > g c (critical value of g) and L ^> £ (the 
correlation length of z), the Binder cumulants reach the values C 2 = 1 and C 4 = 1. 
For < g < g c and L ^> £, we have C 2 = n/2 « 1.57 and C 4 = 3. For < g < g c and 
L <C £, the moments C 2 and C4 vary only weakly with the linear size L. Therefore 
the critical value of g can be estimated from the common intersection of C 2 and C4, 
respectively, as a function of g for several values of L [23], |2"2| . 



3.3 Results and Discussions 

The number of type-a junctions is determined by Eq. (f2.7j) . and we rewrite it as 

where £i is the membrane height at site % at some instance. Thus, we can calculate the 
density of type-a junction in the simulations. When the system reaches equilibrium, 
the dimensionless chemical potential (energy unit is chosen to be ksT) of type-a 
junctions (p, a ), f ree ligands (/Ilo), an d free receptors (jiRa) satisfy the condition 

fia = Jj-RoL + &La- (3-4) 
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The chemical potentials of the free ligands and receptors are related to their densities 
by ipLa = G^ La and ipR a = e^ Ra , and we assume fiRa = j^La- Therefore, the relations 
between -B e //+ and total density of junctions and free ligands and receptors are shown 
in Tab. 13.21 It is clear that the density increase as the effective binding energy 
increases. For higher total effective binding energy ( E e ff+), the total junction density 
is higher. It is because the effective binging energy is the sum of the binding energy 
and the entropy lost due to the formation of ligand-receptor pair (junction). For 
larger effective binging energy, the ligand and receptor prefer to bind to each other. 
In biological systems, typical values of the number of ligands (N La ) and receptors 
(Nlo) are on the order 10 5 per cell, where the area of the cell (A c ) ~ 10 — 10 4 fim 2 . 
In general case, we except that total number of free type-a ligands and receptors 
are the same order as Ni a and Nr u . Thus typical values of ipia and ipR a should 
be O(10 -2 ), and typical binding energy Esa ~ 10-20 ksT, typical E e ff\ and E e ff 2 
should be of order unity, and they can be positive of negative. As shown in Tab. 13.21 
our choice of parameters correspond to typical biological systems. 

Typical C 2 (Fig. 13. 5|) and C4 (Fig. 13 diagrams for several values of L are shown. 
From Fig. 13.51 and 13.61 we find that g c ~ 2.8. We can estimate critical point in this 
way for different E e ff + . Therefore, we plot the critical values of g with respect to 
E e ff + in Fig. 13.71 It shows that the critical value of g decreases with the increase of 
E e ff + , this is because the fluctuation of the membrane is smaller when the junction 
density is higher. 
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Table 3.2: Total density of junctions and free ligands and receptors 



E e ff + junctions(%) free ligands and receptors(%) 

-8 ~ 1.6 ~2x KT 4 

-6 ~ 4.8 ~5x 1CT 3 

-4 ~ 11.2 ~ 10~ 3 

-2 ~ 28 ~ 3 x 1CT 3 

~ 50 ~ 9 x 10~ 3 

2 ~ 73 ~ 0.03 

4 ~ 87 ~ 0.07 

6 ~ 95 ~ 0.2 

8 ~ 98 ~ 0.5 
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Figure 3.5: C 2 as a function of g for L = 16 (square), 32 (circle), and 64 (triangle) 
with parameters A x = A 2 = 0.2, E e ff + = 2 (-E e //i = ^e//2 = 1), and l x = 15. The 
common intersection point in this case is near g = 2.8. 
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Figure 3.6: C 4 as a function of g for L = 16 (square), 32 (circle), and 64 (triangle) 
with parameters A x = A 2 = 0.2, E ef f + = 2 (-E e //i = E ef f 2 = 1), and t\ = 15. The 
common intersection point in this case is near g = 2.8. 
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Figure 3.7: The critical value of g as a function of -E , e //+ f° r ^eff- = 0, A + = 0.4, 
A = 0, and l\ = 15. The critical value of g decreases with the increase of E e ff + , 
this is because the fluctuation of the membrane is smaller when the junction density 
is higher. The data points in the figure are partial result of our simulations. It is 
because the statistical error increases as the value of -E e //+ decreases. Thus, we show 
the data points for —4 < E e ff + < 8. 
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Chapter 4 



Conclusions 

In this thesis, we have studied the adhesion-induced phase separation of multi-component 
membranes by theoretical and numerical methods. In this system, two different types 
of junctions (they are different in the natural length and the flexibility) mediate the 
membrane adhesion. The phase diagram for a system with A > is showen schemat- 
ically in Fig. 14.11 Three methods are used to study the phase separation of this 
system. 

1. In mean field theory, we integrate over all distributions of junctions in the partition 
function and obtain an effective membrane potential. The junction height difference 
and the junction flexibility difference are two important factors that affect the phase 
behavior. We discuss the phase diagrams in both symmetric case (both types of 
junctions have the same flexibility but different natural length) and asymmetric case 
(different type of junctions have different flexibilities and natural length). After an- 
alyzing the properties of the potential, we find that in symmetric case, the phase 
separation is driven by the height difference of the junctions. In the region of suffi- 
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small g region 
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-0.4 -0.3 -0.2 -0.1 
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large g region 



0.1 0.2 0.3 



■eff- 



Figure 4.1: The schematic phase diagram for mean field theory (solid line) and Gaus- 
sian approximation (dashed line). 
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ciently large height difference in the asymmetric case, the junction height difference 
is the main factor that drives the phase separation. However, in the small height 
difference region in the asymmetric case, the minimum energy of the more rigid junc- 
tions is "lowered" by the binding energy of the softer junctions. As a result, in this 
region the phase coexistence occurs when the effective binding energy of the more 
rigid (softer) junctions are smaller (larger). Thus, the phase boundary moves toward 
smaller E e ff_. 

2. Gaussian approximation is used to study the fluctuation effect on the phase 
coexistence curves of the asymmetric case. When the junction height difference is 
large, because the softer junctions can be stretched or compressed easier than the 
more rigid junctions and have higher entropy, comparing to the mean field theory, 
the phase coexistence curves move to higher effective binding energy of the more 
rigid junctions. Furthermore, the softer the junctions are, the higher the the entropy 
they have. Thus, the phase boundary shift is more significant for system with larger 
junction flexibility difference and smaller total junction effective binding energy. In 
small junction height difference region, the phase boundary moves toward the mean- 
field critical point as the junction height difference decreases (because Gaussian theory 
and mean field theory predict the same critical point). 

3. Monte Carlo simulations simulate the symmetric case to study the effect of 
junction density on the critical point. We find that the junction density increases 
when the total effective binding energy increases. Meanwhile, the phase separation 
occurs at larger junction height difference when junction density of the system is 
lower. Because the fluctuation of membrane is smaller when the junction density is 
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higher. 

Experimentally, the transition from shorter-junction-rich phase to longer-j unction- 
rich phase can be induced by applying a mechanical pulling force on the system. This 
is because the effective binding energy of the longer junctions increases and that of 
the shorter junctions decrease under a pulling force. Thus, the phase transition can 
be observed in typical experiments. 

In summary, we have provided a general physical picture for the adhesion-induced 
phase separation of multi-component membranes. In our theory, we assume that 
the membranes bind to each other in the whole process. It is possible to study 
binding/unbinding transition in the future. In this case, the membrane-membrane 
collisions may drive a different kind of phase separation. Another possible future 
work is to study the dynamics of adhesion and detachment of this type of systems. 
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APPENDICES 



Appendix A 



Elasticity of two fluctuating 
membranes 

This appendix shows the calculation that reduces a two-membrane Hamiltonian to a 
single-membrane Hamiltonian. If we consider only the bending energy of the mem- 
branes then the Hamiltonian of the system with two membranes in the solvent can 
be written as 

H = J d2r { T [V ^ l(r)]2 + Y [ y2 ^W] 2 } ' ^ 

where z\ and z 2 are the height of the upper membrane and the lower membrane 
from the reference plane, respectively. Define u = Z\ + z 2 and h = z\ — z 2 , thus 
z x = \{u-h) and z 2 = \{u + h). Substitute [V 2 2i] 2 = [|(V 2 m - V 2 h)] 2 and [V 2 z t ] 2 = 
[|(V 2 m - V 2 h)} 2 into Eq. (jUJ) leads to 

H = l -Jd 2 r[{K l + k 2 ) [V 2 u{r)] 2 + {k x + k 2 ) [V 2 h(v)} 2 
+ 2{k x - K 2 )V 2 u{r)V 2 h{v)} 

42 



2 | 1 KiK 2 

a r 



2 Ki + K 2 



[V 2 Mr)] 2 + ^[V 2 «(r) 



H V h(r) 



(A.2) 



We integrate out the second term in Eq. (jA.2|) and define k = , then the 

Hamiltonian becomes 



H = Jd>r{^[V>h(r)} 2 }. 



(A.3) 



Similarly, for systems with non-vanishing surface tension, 



H 



(A.4) 



where 7 



71 72 
71 +72 ' 
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Appendix B 



Nondimensionalization of the 



Hamiltonian 



The Hamiltonian of the system is 

H = J dxdy [ K - [V 2 h(x, y)T + l [Vh(x, y)f 

2 2 ^ 

+ J2 ^(x, y)^- [h{x, y) - K} 2 - J2 ®*E Ba I . (B.l) 

a=l a=l ) 

The unit length in z-direction is ho = aj ^ jr^f- The discretized Laplacian and gra- 
dient of h in two dimensional space are 

v 2 /i d2fl i &2h 

dx 2 dy 2 

d h x+ a 7 y-h x -a d .h x , y+ «-h Xt y-a 

= t-( — +7T ) 

ox a ay a 

^Q^x+a^y h x ^ a y -\- hx^y+a -\- h x ^y^ a h X y) 



'{^x+a,y ^x—a,y @-x,y—a ^x,y+a ^x,y) 



a v k 



a v k a v k 



' i7 -A, ; /.,, ; -y^-X', (B.2) 
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and 



Vh 



dx dy 



1 

2a 
1 



{hx+a,y h x —a,y)l \hx,y+a h Xl y—a)^ 



2a V k 



[(^x+a,y @-x-a,y)~X- + (£x,y+a ^x,y-a)y] 



1 '^ V4 ,„ = ^,/^ w ,. 



2a V /c 



2a V k 



(B.3) 



Also defined here is <f> a (x, y) = a 2 $ a (r), the number of type-a junctions at site (x, y). 
Hence, the non-dimensionalized Hamiltonian of the system is 



i=i ^ 

2 



Q = l 



1 X ^7 2 



Bo 



(B.4) 
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